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ABSTRACT 

The Palomar Transient Factory (PTE) provides multiple epoch imaging for a large fraction of the celestial 
sphere. Here we describe the photometric calibration of the PTE data products that relates the PTE magnitudes 
to other magnitude systems. The calibration process utilizes Sloan Digital Sky Survey (SDSS) r ^ 16 mag 
point source objects as photometric standards. During photometric conditions, this allows us to solve for the 
extinction coefficients and color terms, and to estimate the camera illumination correction. This also enables 
the calibration of fields that are outside the SDSS footprint. We test the precision and repeatability of the PTE 
photometric calibration. Given that PTE is observing in a single filter each night, we define a PTE calibrated 
magnitude system for /?-band and g-band. We show that, in this system, « 59% (47%) of the photometrically 
calibrated PTE /?-band (^-band) data achieve a photometric precision of 0.02-0.04 mag, and have color terms 
and extinction coefficients that are close to their average values. Given the objects' color, the PTE magnitude 
system can be converted to other systems. Moreover, a night-by-night comparison of the calibrated magnitudes 
of individual stars observed on multiple nights shows that they are consistent to a level of « 0.02 mag. Most 
of the data that were taken under non-photometric conditions can be calibrated relative to other epochs of the 
same sky footprint obtained during photometric conditions. We provide a concise guide describing the use 
of the PTE photometric calibration data products, as well as the transformations between the PTE magnitude 
system and the SDSS and Johnson-Cousins systems. 
Subject headings: techniques: photometric - catalogs 



L INTRODUCTION 

The Palomar Transient EactorvPI (PTE; Law et al. 2009; 
Rau et al. 2009) is a synoptic survey designed to explore 
the transient sky. The project utilizes the 48" Samuel Oschin 
Schmidt Telescope on Mount Palomar. The telescope has a 
digital camera equipped with 11 active CCDs, each 2Kx4K 
pixels (Rahmer et al. 2008). Each PTE image covers 
7.26 deg^ with a pixel scale of 1.01"pix^^ Erom the begin- 
ning of the PTE survey in March 2009 until January 2011, 
most of the images were taken with the Mould R filter. Start- 
ing January 2011, we performed the PTE main survey in ^- 
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band during dark time and in /?-band during bright time. In 
addition, a few nights around times of full Moon are used for 
taking images with narrow-band Ha filters. A PTE system 
overview and review of the first year's performance are given 
in Law et al. (2010). 

Accurate photometric calibration is a non-trivial task since 
one must know both the atmospheric transmission (e.g., Pad- 
manabhan et al. 2008; Burke et al. 2010) and the op- 
tical/detector system response as a function of wavelength, 
time, and sky position. Moreover, observations of flux stan- 
dards are required, and the true spectral energy distribution of 
these standards needs to be known. In practice simplifying 
assumptions are made in order to achieve solutions which are 
good to the few percent level. Eor example, in many cases, it 
is customary to use the calibrated magnitudes of the standard 
stars instead of their spectral energy distribution, and to as- 
sume that the atmospheric transmission is a smooth function 
of airmass. We note that relative-photometry approaches (e.g., 
GilHland & Brown 1988; GilHland et al. 1991; Honeycutt 
1992; Everett & Howefl 2001; Ofek et al. 201 la) can achieve 
better accuracy from the ground, but are limited by scintilla- 
tion noise (e.g.. Young 1967; GilHland & Brown 1988), Pois- 
son noise and flat-fielding errors. 

There are several approaches to photometric calibration. 
Eor example, the Sloan Digital Sky Survey (SDSS; York et 
al. 2000) images are calibrated using an auxiliary 20-inch 
telescope, which determines the photometric condition on a 
nightly basis (Hogg et al. 2001). This telescope measures 
the extinction and photometric zero point using a network of 
standard stars (Smith et al. 2002), which, in turn, are tied to 
the standard star BD+17:4708. Furthermore, the photometric 
uniformity in the SDSS is achieved by the algorithm described 
in Padmanabhan et al. (2008). The SDSS photometry is uni- 
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form to better than 2% in all bands (Adelman-McCarthy et al. 
2008). 

Ofek (2008; see also Pickles & Depagne 2010) has sug- 
gested using Tycho-2 (H0g et al. 2000) stars to photometri- 
cally calibrate astronomical images. The SDSS griz magni- 
tudes of such stars are based on synthetic magnitudes of stel- 
lar spectral templates fitted with the Hipparchos Bj, Vt mag- 
nitudes and the 2MASS (Skrutskie et al. 2006) JHK mag- 
nitudes. However, this strategy requires images containing 
unsaturated Tycho calibration stars which lie in the w 9-12 
magnitude range. The saturation limit for the PTF g- and R- 
band survey is typicall}(3 around 14 magnitude. 

Herein, we describe the method we have developed for the 
photometric calibration of PTF data taken in the g and R 
bands. The photometric calibration of the PTF Ha survey will 
be described elsewhere. We note that the relative photometric 
calibration of PTF currentl\F^ achieves precision as good as 
3mmag, in given fields, at magnitude 15 (e.g., van Eyken et 
al. 2011; Agiieros et al. 201 1; Law et al. 201 1; Levitan et al. 
201 1; Polishook et al. 2011). The PTF relative photometry 
pipeline will be described in Levitan et al. (in prep.). 

Our method of photometric calibration for PTF data is sim- 
ilar to the "classical" method of observing a small number 
of calibration stars (e.g., Landolt 1992) through various air- 
masses and assuming photometric conditions - i.e., the atmo- 
sphere transmission properties are constant in time and are a 
continuous function of airmass. The main difference is that 
we are using SDSS stars as standard stars and we typically 
observe 10^ SDSS stars with high S/N in each CCD per 
night. Another important difference is that PTF observations 
are done in a single filter each night. Therefore, in order to 
relate the PTF calibrated magnitudes to one of the common 
absolute systems one need to apply color terms. The paper is 
organized as follows. In 5j2] we describe the PTF photometric- 
calibration method. In !j3] we discuss the illumination correc- 
tion, while Sj] describes a related problem that plagues early 
PTF data. In Sj5]we describe the data products and what cal- 
ibration information is stored with the image-product data. 
The performance of the PTF magnitude calibration is given in 
Sj6] and the derived photometric parameter statistics are dis- 
cussed in SjT] We provide the color transformation between 
the PTF magnitude system and other systems in [Js] Finally, 
we conclude the paper in ij9] Unless specified otherwise, the 
statistics given here are based on all PTF data obtained from 
Mai-ch 2009 to July 2011. 

2. PHOTOMETRIC-CALIBRATION METHOD 

PTF has two main data-reduction pipelines. The first is 
for real-time (ss 30min; e.g.. Gal- Yam et al. 2011) image 
subtraction and transient detection, hosted by the Lawrence 
Berkeley National Laboratory (Nugent et al., in prep.). This 
paper will make no further mention of this pipeline, as only 
the second pipeline is relevant here, as discussed below. 

The photometric calibration described in this paper is im- 
plemented in the second pipeline. This pipeline, hosted by the 
Infrared Processing and Analysis Center (IPAC), performs fi- 
nal image reduction and extracts the source catalogs. The pro- 
cessing includes splitting the multi-extension FITS images. 

The PTF camera electronics were modified on 2009 Oct 22 to increase 
the dynamic range. Before this date, the saturation limit was around 15th 
magnitude. 

For relative photometry, we use a scheme similar to that proposed by 
Honeycutt (1992) with some modifications outlined in Ofek et al. (201 la). 



de-biasing, flat fielding, astrometric calibration, generation of 
mask images, source extraction and photometric calibration 
(the subject of this paper). This pipeline is described in Grill- 
mair et al. (2010) and Laher et al. (in prep.). 

Our photometric-calibration process runs on PTF data sepa- 
rately for each night, filter and CCD. It first attempts to match 
the sources extracted from the PTF images taken during a 
given night with SDSS-DR7 PhotoPrimarjH point sources 
(i.e., SDSS type= 6). In order to assure good photomet- 
ric quality, only SDSS stellar objects with photometric errors 
smaller than 0.05 mag in r- and /-bands and which are fainter 
than 15 mag (to avoid saturated stars) in the g-, r- and /-bands 
are used. The photometric solutions are calculated only if 
more than 30 science images were taken during the night, and 
only if we were able to select more than 1000 SDSS stars in 
all the images taken with a given CCD and filter during the 
entire night. 

We use the SDSS-matched stars as a set of standard stars 
and solve for the photometric zero points, airmass terms and 
color terms in a given night. The fitting process is done sep- 
arately for each one of the 11 active CCDs and the g- and 
/?-band filters. This is required since the CCDs are not iden- 
tical and some of them have a different spectral response (see 
Law et al. 2009 for details). For observations taken using the 
7;-ban43 filter we fit the following model: 

''SDSS - ^PTF = + '^c,^ ('"SDSS - 'SDSS ) 

+ aa.RAM+ aac,MM(rsDSS - i'SDSs) 
+ CCt,R{t - tra) + OaAt - tmf- 

-2.51ogio(5f), (1) 
while for g-band observations we fit: 

gSDSS - gpTF = '^Pg + «c.g fesDSS - /"SDSS ) 

+a^^gAM + aac,gAM(gsDSS - '"sdss) 

-2.51ogio(5f). (2) 

Here, /p" p is the PTF instrumental magnitudes in band / (ei- 
ther R or g), fsDSS is the SDSS magnitude in band / (either 
g, r or /), ZPf is the photometric zero point for filter /, a^j 
is the color term for filter /, tta./ is the extinction coefficient 
(airmass term) for filter /, Ctac./ is the airmass-color term for 
filter /, AM is the airmass, at./ and Oaj are the polynomial 
coefficients for the change in the zero point of filter / as a 
function of time, t in days, during the night, where in days 
is the middle of the night, and 5f is the exposure time in sec- 
onds. 

The instrumental magnitudes used in the photometric cali- 
bration process are based on the SExtractor (Bertin & Arnouts 
1996) MAG_AUTO magnitud^3 (see also ^O ), with an in- 
ternal SExtractor zero point of 0. The above set of equations 
are solved using linear least-squares fitting. The errors in 
'■SDSS-Z^ii'^;. are taken as (Argoss^ + A/?j?^^^ + 0.015^) 
where ArsDSS is the SDSS magnitude error, ARJ^^p is the PTF 
magnitude errors and 0.015 is the assumed internal accuracy 
of the SDSS photometric calibration (Adelman-McCarthy et 
al. 2008). The fit is performed iteratively, up to 3 times, with 

PhotoPrimary is an SDSS table. See definitions in the SDSS Schema 
Browser: http://cas.sdss.org/dr7/en/help/browser/browser.asp 

The PTF Mould R filter is similar in shape to the SDSS r-band filter, but 
shifted 27 A redward. 

" Defined with kron_fact =1.5 and min_radius = 2.5. 
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sigma clipping of 3(7. The sigma clipping ensures removal 
of stars with bad photometry (e.g., influenced by cosmic rays; 
saturated pixels). 

Next, in order to be able to correct for zero-point variations 
across a given CCD (see the residuals from the best fit of 
Equations [T] or |2] are binned in cells of 256 x 256 pix^ along 
the X and Y dimensions of each CCD. In each cell, we take 
the mean of the residuals and subtract from it the mean of the 
residuals in the center-of-image celQ in order to render the 
residuals relative to the residual at the image center The re- 
sulting coarse image of the mean of the residuals is linearly 
interpolated to generate an image at the resolution of PTF im- 
ages of pixel-to-pixel zero point variations. This image prod- 
uct is equivalent to an illumination correction, and herein also 
called the zero point variation map (ZPVM; see f|3j. The un- 
certainty in the ZPVM is estimated by calculating the standard 
deviation (StD) of the residuals in each celFT. 

We note that the PTF illumination correction images are 
usually smooth on large scales (comparable with the CCD im- 
age size). Therefore, in most cases the ZPVM can be repre- 
sented by low-order polynomials. In order to provide users 
with a simpler version of the illumination correction, we also 
fit versions of Equations [T] and |2] that consists of a low-order- 
polynomial representation of the ZPVM, e.g., for the 7?-band: 

''SDSS — ^PTF = + ^c.R [rSDSS — 'SDSS ) 

+aaMM + a^cRAMirsDSS - 'sdss) 

+ at,R{t - fm) + at2.R{t - tmY 

-2.51ogio(5f) 

+ 0!xi.ff (X — 2^size)/^size 
+ ~ 2^size)/Jsize 

+ %2,/f[(J' ~ 5J'size)/isize]^ 
+ %3,«[(J' — 2^size)/Jsize]^ 
+ OJxy,«(^ ~ 2^&i7e){y ~ 2^size)/ (''^size^size)(3) 

where X and Y are the positions on the CCD, Xsize is the size 
of the CCD, in pixels, in the X dimension (2048 pix) and Fsize 
corresponds to the Y dimension (4096 pix). A similar equa- 
tion is fitted for the g-band observations (see Equation |2] for 
analogy). Here, the illumination correction is represented by 
polynomials with coefficients ttxi./, 0Ly2.f, CtyS./, and 

ttxy./- This representation usually provides a good estimate of 
the illumination correction over the image. We note that for 
each solution we store a variety of information regarding the 
quahty of the fit (see [jSj. 

3. ILLUMINATION CORRECTION 

The ZPVM and its polynomial representation amount to 
an illumination correction that describes the variations in the 
photometric zero point spatially across the image (e.g., Reg- 
nault et al. 2009). Naively, such zero-point variations should 
be removed by the flat-fielding process. However, for exam- 
ple, if the fraction of flux in the stars' point spread function 
(PSF) wings, relative to the flux integrated to infinity in the 
PSF, is variable as a function of position on the CCD then this 
may induce variations that are not removed by the flat-fielding 
process. Such CCD position-dependent variations are indeed 
detected in PTF images (e.g., fj4]l. Figure [T] presents two 
examples of illumination correction images. The left panel 

The position of the center of thi.s central cell is x = 1025, v = 2049. 
^' The StD is not divided by the square root of number of data points in 
each cell. Therefore, it represents the scatter rather than the error in the mean. 
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Fig. 1 . — Two examples of ZPVM images. The color coding represents the 
vaiiations in photometric zero point with respect to the center of image. Left 
panel: a good case with peak-to-peak range of 0.016 mag in the zero point 
vaiiations across the image. This ZPVM image was constructed based on 
data obtained on 2011 July 13 with CCDID= 8 and the R-band filter. Right 
panel: One of the worst examples with peak-to-peak range of 0.22 mag in 
the zero point across the image. This ZPVM image was constructed based on 
data obtained on 2010 May 7 with CCDID= 5 and the R-band filter. In such 
bad cases the polynomial representation may deliver poor quality results. We 
note that this bad example is not representative of PTF data and that typical 
images obtained using this system have good image quality. 

shows a typical case, in which the illumination corrections 
have a low range variations. The right panel shows one of the 
worst cases we encountered thus far, with larger amplitude 
variations (see S|4|. We note that the median peak-to-peak 
variations in the zero point across the image over all photo- 
metric nights, filters and CCDs is about 0.06 mag. 

In the first year of the PTF project following first light, we 
detected a problem with the PTF image quality. This issue is 
referred to as the fogging problem, and it manifested itself as 
a diffuse halo around bright stars and a prominent large scale 
non-uniformity in the illumination corrections. The fogging 
problem and our testing and subsequent apparatus modifica- 
tions for ameliorating its effects are described in S|4] 

4. THE FOGGING PROBLEM 

The PTF camera window, which also acts as a field flattener, 
is fairly large (« 320 cm^). The window surface temperature 
is w 16 °C below the ambient temperature. To prevent con- 
densation from forming on the window, a constant stream of 
dry air or nitrogen is blown across the window's outer sur- 
face. In the original installation of the camera, new nylon 
tubing was used to convey the gas to the window. This tub- 
ing was selected for its reported high durability, flexibility and 
low outgassing properties. Once the window fogging problem 
was discovered, a test chamber was assembled, which sim- 
ulated the cold window environment in which dry nitrogen 
could be conveyed across a cold test surface. 

The outgassing properties of various tubing materials were 
evaluated in the test chamber. The nylon tubing used in the 
original camera installation was found to have a moderate 
amount of outgassing which produced volatiles that would 
condense on the cold test surface. Tubing made from Teflon 
FEP (fluorinated ethylene propylene) was found to have low 
outgassing with no detectable volatile condensation on the 
cold test surface. The camera window was cleaned and all 
tubing between the dry nitrogen supply and the camera win- 
dow was replaced with this material on 2010 September 2. 
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TABLE 1 

PTF CAMERA WINDOW CLEANING DATES 



Year 


Cleaning date 
Month 


Day 


Comments 


2009 


07 


17 




2009 


07 


29 




2009 


08 


11 




2009 


08 


20 




2009 


09 


15 




20U9 


uy 


29 




2009 


10 


21 




2009 


11 


10 




2010 


01 


11 




zUlU 


01 


L 1 




2010 


02 


04 




2010 


02 


18 




2010 


03 


04 




2010 


04 


06 




2010 


04 


21 




2010 


05 


12 




2010 


06 


02 




2010 


06 


17 




2010 


07 


01 




2010 


08 


04 




2010 


08 


16 




2010 


09 


02 


new tubing 



Note. — Local dates the camera window was cleaned (e.g., the last cleaning date is 
2010 Sep 02, so all the data taken on UTC 2010 Sep 03 and afterward is affected by this 
cleaning). At these dates, the illumination correction may change markedly. 



No window fogging issues have been observed since the tub- 
ing material was replaced. 

We note that the camera window was cleaned several times 
during PTF operations, and after each cleaning an immediate 
reduction in ZPVM deviations from zero is expected. The 
cleaning dates are listed in Table [T] Fortunately, the 

method described in |2]provides an estimate of the zero-point 
variations across the CCD image. Therefore, the illumination 
correction is used to correct our photometry, and to monitor 
the PTF image quality. 

5. CALIBRATION PRODUCTS 

The parameters required for calibrating the instrumental 
magnitudes are stored in the PTF database, as well as in the 
image and catalog FITS file headers (! 15.1I ). In addition, a full 
version of each illumination-correction image is stored as an- 
cillary calibration data in a FITS file, and associated metadata 
are stored in the PTF database (i )5.2l i. The use of these data 
products is discussed in i]5.3l 

5.1. FITS headers 

The parameters based on the solution of Equation [3] are 
stored in the image headers. Moreover, generated from each 
image is a FITS binary table containing all the sources ex- 
tracted from that image using SExtractor The FITS binary 
table catalogs contain three header + data units (HDUs) and 
two FITS extensions. The first HDU (i.e., the primary HDU) 
contains a header with metadata from SExtractor and no data. 
The second HDU, or first FITS extension, is a binary table 
containing the properties of the sources extracted from the 
processed image, with its header giving information about 
the table columns. The source properties in the binary ta- 
ble include aperture-photometry instrumental magnitudes and 
other SExtractor outputs. The binary table has also been aug- 
mented with an additional column called ZEROPOINT, which 
stores the source-dependent zero points computed from our 



photometric-calibration parameters; these zero points already 
include the 2.5 log(5f ) contribution for normalizing the image 
data by the exposures time, 5t, in seconds, and so straightfor- 
wardly adding the instrumental magnitudes to their respective 
zero points will result in calibrated magnitudes. The third 
HDU, or second FITS extension, contains an empty image 
whose header is a copy of the processed-image header, which 
includes information about the photometric calibration as pre- 
viously mentioned. The header keywords related to the pho- 
tometric solutions are summarized in Table |2] 

The FITS headers contain additional parameters related to 
photometric calibration. However, these parameters are based 
on different routines which are no longer supported. For com- 
pleteness, these parameters are listed in Table|2]under the ob- 
solete parameters section. These keywords are products of our 
first attempt at photometric calibration of PTF images, which 
is less sophisticated than the method presented here. Both 
methods are still being run in the pipeline and the results from 
both methods are still being written to the FITS headers of the 
processed images. A detailed description of the first version 
of photometric calibration is beyond the scope of the current 
paper, but will be available from online documentation. 

5.2. The illumination correction FITS image 

For each night, CCD and filter, we generate two additional 
calibration products, and package them in separate single- 
extension FITS -image files: (i) the illumination correction 
image (or ZPVM image); and (ii) an image of the standard 
deviation of the illumination correction. Both images are cre- 
ated by linearly interpolating the values of the means and stan- 
dard deviations of the residuals from Equation[T]or|2](see fj2]i. 
The ZPVM image and associated standard deviation image 
are calculated in units of magnitude. The values of the ZPVM 
image are shifted by a constant so that they are equal to zero 
in the image center Setting the value at the image center to 
zero makes the solutions based on Equations [T] or [2] similar to 
those based on Equation[3] 

The illumination correction images are not applied to the 
individual processed PTF images. One reason is that the illu- 
mination correction images are reliable only for photometric 
nights. The selection of photometric nights is discussed in ^ 

5.3. Use of data products 

Each processed PTF image is accompanied by an associ- 
ated FITS binary table containing a catalog of objects ex- 
tracted from the processed image by SExtractor. The catalog 
source magnitudes are instrumental magnitudes, and users of 
this product can convert instrumental magnitudes into cali- 
brated magnitudes using the provided PTF photometric pa- 
rameters and the following procedure. In order to convert the 
MAG_AUTCE3 instrumental magnitudes to calibrated magni- 
tudes, one needs to apply Equation[T]or|2]and add the value of 
the ZPVM image at the location of the object. Alternatively, 
it is possible to directly use Equation[3] 

Applying Equation[Tl|2]or[3]requires the colors of individ- 
ual sources in the SDSS system. If the color of a source is 
known, from PTF observations or other surveys (e.g., SDSS), 
then it is possible to use these equations to get the magni- 
tude of that source in the SDSS system. However, not all PTE 
fields were observed in both the g and R filters or are within 

For description of the MAG_AUTO magnitude see the SExtractor man- 
ual : http ://w w w. astromatic .net/ soft ware/sex tractor 
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TABLE 2 

The photometric calibration FITS header keywords 



Keyword Symbol Units Description 



FILTER 


/ 






Filter name (e.g., g, R) 






AFXPTTMF 


St 






Actual exposure time 






APSFILT 








SDSS filter used in photometric calibration (e.g. 


g, r) 


APSCOL 








SDSS color used in photometric calibration (e.g 




— r, r — (') 


APRMS 




mag 




RMS of the residuals in photometric calibration 






APBSRMS 




mag 




RMS of the residuals in photometric calibration for stars brighter than 16th magnitude in r-band. 


APNSTDIl 








Number of calibration stai's in first iteration 






APNSTDIF 








Number of calibration stai's in final iteration 






APCHI2 








7 ' of photometric calibration 






APDOF 








Number of degrees of freedom in photometric calibration 


APMEDJD 


t,„ 


day 




Median JD for the night 






APPNOl 








Name of parameter photometric calibration 01 (i 


.e., 


'ZeroPoint') 


APPAROl 


ZPf 


magADU"' s"' 




Value of parameter photometric calibration 01 






APPARbOl 


tsZPf 


mag 




Error in parameter photometric calibration 01 






ArriNUz 






Name of parameter photometric calibration 02 (i 


.e.. 


'ColorTerm') 


A I ) 1 1 A Tj no 
ArrAKUZ 


"c/ 


mag mag^ ' 




Value of parameter photometric calibration 02 






A Tin A m~7m 

APPARbUz 


AOc./ 


mag mag^ ' 




Error in parameter photometric calibration 02 






APPJNUj 






Name of parameter photometric calibration 03 (1 


.e.. 


'AirMassTerm') 


A DDA T?ni 




mag airmass^' 




Value of parameter photometric calibration 03 






APPARbOj 


Atta,/- 


mag airmass^' 




Error in pai'ameter photometric calibration 03 






APPN04 






Name of parameter photometric caUbration 04 (1 


e.. 


' AirMassColorTerm ' ) 


A DDA 'DCiA 


«ac,/ 


mag mag airmass 


Value of parameter photometric calibration 04 






APPARE04 


Attac./ 


magmag^' airmass"' 


Error in parameter photometric calibration 04 






APPJNU5 








Name of parameter photometric calibration 05 (i 


.e.. 


'TimeTerm') 


A T3PA Dn^ 




mag day ' 




Value of parameter photometric calibration 05 






A nnA TiTTrtc 

APPARh05 


AOt./ 


mag day"' 




En'or in parameter photometric calibration 05 






APPJNUo 






Name of parameter photometric calibration 06 (i 


.e.. 


'TimeTerm2') 


A PPA Dn^ 




mag day"^ 




Value of parameter photometric calibration 06 






ArrAREUO 


Aa,2./ 


mag day"^ 




Error in parameter photometric calibration 06 






A PPMOT 

ArrJNU / 






Name of parameter photometric calibration 07 (i 


.e.. 


'XTerm') 


A PPA pn7 
Ar rAKU / 


«xL/ 


mag(Xs,ze pix)"' 




Value of parameter photometric calibration 07 






A PPA pr;n7 
Ar rAKllU / 


Aoixi,/ 


mag(Xsize pix)"' 




Error in parameter photometric calibration 07 






A ppwnfi 
Ar rlNUo 






Name of parameter photometric calibration 08 (1 


e.. 


'YTerm') 


APPAR08 


«yl./ 


mag(y,i2epix)"' 




Value of parameter photometric cahbration 08 






APPARE08 


Aayi.y 


mag(ysizepix)"' 




En'or in pai'ameter photometric calibration 08 






APPN09 






Name of parameter photometric calibration 09 (i 


.e.. 


'Y2Term') 


APPAR09 


072./ 


mag (Fsije pix)"^ 




Value of parameter photometric calibration 09 






APPARE09 


AOyO f 


mag (y,i2e pix)"2 




Error in parameter photometric calibration 09 






APPNIO 






Name of parameter photometric calibration 10 (i 


.e.. 


'Y3Term') 


APPARIO 




mag (y,i2e pix)"3 




Value of parameter photometric calibration 10 






APPAREIO 


Aay3 ^■ 


mag (Fsize pix)"^ 




Error in parameter photometric calibration 10 






APPNll 






Name of parameter photometric calibration 1 1 (i 


.e., 


'XYTerm') 


APPARll 


"xy./ 


mag (Xs,zeFsize)" 


. -2 
pix 


Value of parameter photometric calibration 1 1 






APPAREll 


AOxy,/- 


mag(Xsi2e}'sizc)"' 


. -2 
pix 


Error in parameter photometric calibration 1 1 






obsolete parameters 














PHTCALEX 








Flag indicating if old photometric calibration routines were executed 


PHTCALFL 








Flag for image is photometric (0=N, 1=Y) 






PCALRMSE 




mag 




RMS of the residuals from (zeropoint, extinction) data fit. 


IMAGEZPT 




mag 




Image magnitude zeropoint 






COLORTRM 




mag 




Image (e — r) color term 






ZPTSIGMA 




mag 




Robust dispersion of SExtractor— SDSS magnitudes 


IZPORIG 








Photometric-calibration catalog origin (e.g., SDSS) 




ZPRULE 








Photometric-calibration method (e.g., DIRECT) 






MAGZPT 




mag 




Magnitude zeropoint at airmass= 1 






EXTINCT 




mag airmass"' 




Extinction 







Note. — FITS header keywords repre,senting the nightly best-fit photometric calibration. In all cases the parameters are based on fitting Equationl3l Symbol represents the parameter 
designation used throughout this paper (e.g., Eqs0[3}. The set of keywords below the horizontal line are obsolete parameters (see text). 



the SDSS footprint. In these cases, setting the color terms 
and Uac.f to zero will return the object's magnitude in the 
non-color-corrected PTF magnitude system (i.e., the natural 
PTF magnitude system). We denote magnitudes in this sys- 
tem by /pTF> where / is the filter moniker. PTF magnitudes 
which are color-corrected (i.e., the source color is known and 
applied to Equations [Tl|2] or O are denoted by /ptf/sdss- We 
note that applying Equations [T] |2] or [3] is done by isolating 
/sDSS and replacing it with /ptf/sdss in cases where the color 



is known, or /ptf if the colors are set to zero. We note that 
the calibrated /ptf magnitudes of all sources detected in PTF 
images are available from the PTF database. 

As in the case of any other magnitude system, in order to 
convert the PTF non-color corrected magnitudes to other sys- 
tems, one needs to know the color index of the source. In 
some cases, it is possible to use some prior knowledge about 
the mean color index of a class of objects (e.g., RR Lyr stars; 
asteroids), in order to convert their PTF magnitudes (in a sin- 
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Fig. 2. — rsDSS ^^PTF/SDSS 3s 3 function of SDSS r-band magnitude 
(rsDSs) for point sources in a single PTF CCD image. The image was taken 



on 2010 March 18.3696 with CCDID= 
on PTF field 2961. 



were the PTF camera was centered 



gle band) into SDSS magnitudes. For example, 68% of main- 
belt asteroids have r — i color index in a narrow range of about 
0.2 mag (e.g., Ivezic et al. 2001). Therefore, given the PTF 
R-hmd color coefficients, oCc.r ~ 0.2, the scatter in the mean 
color correction for asteroids is about 0.04 mag. Neglecting 
the airmass-color terms (which are typically small; see Ta- 
bleO, the relations between the PTF and the SDSS magnitude 
systems, for objects with stellar-like spectra, are given by: 



and 



''SDSS -^PTF — OCc,R{rSDSS - 'SDSs) 



gSDSS -^FTF = OCc.g(gSDSS " '"SDSs), 



(4) 



(5) 



where the values of the color terms are listed in Table [3] 

In order to apply the illumination correction to individual 
images, we need to multiply, pixel by pixel, the processed 
image by jq^'^^^^^^'*', where ZPVM is the illumination cor- 
rection image. Alternatively, one can use the less accurate 
polynomial representation of the ZPVM image. 

6. PERFORMANCE 

The accuracy of our calibration process, relative to SDSS, 
is discussed in il6.1l while in we discuss the repeatability 
of our photometric calibration. 

6.1. Accuracy 

As described in the PTF calibrated magnitudes can be 
given in two systems. A color-corrected magnitude system, 
denoted by /^ptf/sdss or ^ptf/sdsS' which the known color 
of the object is applied to Equations [T]-|2] or [3] By construc- 
tion, this magnitude system is similar to the SDSS system. 
The other possibility is a "natural system" in which all the 
color terms involving a^.f and ttac./ are set to zero (i.e., the 
color indices of the objects are ignored). These magnitudes 
are denoted by Rpjp or gpjp. 

To demonstrate the current accuracy of the PTF photomet- 
ric calibration, we show in Figure |2] the differences between 
SDSS magnitudes and color-corrected PTF /?-band magni- 
tudes (^pTF/SDSs)' as a function of SDSS magnitude. In 
this particular example, the StD in the rsoss ^ RpTF /sdss for 
bright stars is smaller than k, 25mmag. Figure |3] shows the 
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Fig. 3. — Same as Figure|2]but without applying the color terms to the PTF 
magnitudes (i.e., assuming all the stars have r — i = mag). 
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Fig. 4. — The g — r and r — i color distributions of all SDSS-DR7 point 
sources with r-band magnitude between 16 and 17. 



same, but for PTF magnitudes not corrected for the color term 
(Rpjp). The scatter in the uncorrected color plot is larger, as 
expected. The larger scatter in Fig.[3]is due to the fact that dif- 
ferent stars in the field of view have different colors. The 95th 
percentile range of r — i color of SDSS stars in the r-band 
magnitude range of 16 to 17 is between —0.2 and 1.2 mag 
(see Fig. m. Given that the PTF mean r — i color term is about 
0.2 mag, this corresponds to a color correction between —0.04 
and 0.24 mag. 

We note that in Figures [2l[3] the magnitude difference 
is systematically tending toward larger negative values, 
starting around magnitude 19.5. This is due to the use 
of MAG_AUTO. Indeed, using aperture magnitude (i.e., 
MAG_APER), we do not see any evidence for this bias. 
This is because MAG_AUTO adjusts the aperture used to 
extract the source magnitude for each object. Therefore, 
if MAG_AUTO is used for sources near the survey detec- 
tion limit it is recommended to estimate the aperture cor- 
rection for these sources in order to obtain their unbiased 
magnitudes. We note that the MAG_AUTO is different from 
aperture or PSF magnitudes. In order to convert between 
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MAG_AUTO and other kinds of magnitude estimators, one 
needs to find the aperture correction between these magni- 
tudes. Both MAG_AUTO and aperture magnitudes in prede- 
fined five apertures are available in the FITS binary source 
catalog associated with each image. Future versions of the 
PTF photometric pipeline will use aperture or PSF magni- 
tudes. 

6.2. Repeatability 

The stability over time of the PTF /?-band photometric so- 
lutions is demonstrated in Figure |5] This figure shows the 
scatter in the calibrated magnitude of each star, over multiple 
epochs, as a function of magnitude. The scatter is calculated 
using the 68th percentile range divided by two (i.e., equivalent 
to StD in case of a Gaussian distribution). This is based on 27 
images of PTF fielcS 100037 €00110= taken between 
March and October 2010, which have photometric calibration 
bright-star root-mean squares (RMS) value (parameter APB- 
SRMS in Table|2ll of less than 0.04 mag. 

We choose to use the 68th percentile range instead of StD 
since it is more robust to outliers. For example, if a night was 
"photometric" for 90% of the time (e.g., clouds entered only 
toward the end of the night), then our pipeline might claim 
that the night was photometric, but the calibration of some of 
the data will be poor. Therefore, in order to get the calibrated 
source magnitudes, it is important to average the data taken 
over several photometric nights. We note that we are using 
this approach for the compilation of the PTF photometric cat- 
alog (Ofek et al., in prep.). In turn, this catalog will be used 
in later phases to calibrate the entire PTF data. 

Figure |5] suggests that, at the bright end, the repeatability 
of PTF calibrated magnitudes is good to a level of 0.02 mag. 
Figure |6] shows the same, but for stars in 100 representative 
PTF fields and all CCDs. The typical repeatability of PTF 
magnitudes is field-dependent. The median repeatability over 
all these 100 fields, for all CCDs, for stars between 15 to 
15.5 mag, is about llmmag and the 95-percentile range is 
between 2mmag to 44mmag. We note that PTF g-band ob- 
servations achieve similar repeatability. 

7. PHOTOMETRIC PARAMETER STATISTICS 

The term "photometric night" is not well defined and de- 
pends on the required accuracy. Figure [7] shows the cumula- 
tive histogram of the RMS of the best-fit residuals of bright 
stars (parameter APBSRMS in Table |2| in all the images 
taken by PTF. This figure suggests that about 62% (50%) 
of the /?-band (^-band) images taken by PTF have photo- 
metric calibration with bright-star RMS values smaller than 
about 0.04 magnitude. Figure [8] shows the /?-band airmass 
term (a^.R) vs. the color term {Oc.r), while Figure |9] shows 
the same for g-band. These figures suggest that some 
nights that have small bright-star RMS values (i.e., smaller 
than 0.04 mag) have anomalous airmass or/and color coeffi- 
cients. Therefore, it is better to select data taken in nights for 
which the airmass and color terms are close to their average 
values (e.g., see TableO. Here we define photometric nights 
as nights with bright-star RMS value smaller than 0.04 mag 
and Ucj within the range (Oc./) ±3A[ac,/] and o^r within 
the range (a^j) ±3A[aa,f]- Here, {Oc.f), A[ac.f], {(X^j), and 

23 pjp gg[jj (denoted by PTFFIELD) is a an internal index associated witli 
a unique field position on the celestial sphere, where PTF fields are mostly 
non-overlapping. 

CCD number is denoted by CCD ID and runs from to 11. 
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Fig. 5. — R-band magnitude vs. scatter for stars in PTF field 100037 
CCDID= 0. The plot shows the 68th percentile range divided by 2 as a 
function of the PTF mean magnitude (not color corrected). The 68th per- 
centile range, for each source, is calculated over 27 images of this field taken 
between March and October 2010 which have photometric calibration bright- 
star RMS value (see definition for parameter APBSRMS in Table|2) of less 
than 0.04 mag. We note that the increase in scatter at the bright end is due to 
saturated stars. 
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<Rp-^P> [mag] 

Fig. 6. — Like Figure|5] but for stars in 100 representative PTF fields/CCDs 
which were observed at least on three photometric nights. 

A[aa,/] are the median and \-a width of the color and extinc- 
tion coefficients for each CCD as listed in Table |3] We find 
that 59% (47%) of the R-haxvA (g-band) PTF data were taken 
in such photometric nights. 

The atmospheric extinction coefficients depend on the at- 
mospheric conditions. Therefore, they may vary with season. 
Figure[TO]shows the median 7? -band airmass coefficient (tta.s) 
and the color term {oLcr) as a function of month in the year 
The error bars are given by the StD within each month di- 
vided by the square root of the number of photometric nights 
in each month. This figure indicates that there are small, but 
significant, seasonal variations in both the 7?-band extinction 
and color terms. 

Table [3] presents the median and 68th percentile range 
divided by 2, of the color-term, extinction-coefficient and 
airmass-color terms for each CCD and filter. These are cal- 
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Fig. 7. — Cumulative histogram of tlie bright-star RMS parameter (i.e., 
observations that have bright-star RMS scatter below a particular value) for 
all PTF R-band (red) and g-band (blue) images. 
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Fig. 9. — The same as Figure[8]but for ^-band. The bimodality is due to 
different blue-side spectral response of the various CCDs in the PTF camera 
(see Table[3). We note that most of the PTF data were taken in R-band. 
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Fig. 8. — Airmass term vs. color term of all the PTF R-band data in 
all CCDs. Black points represent images taken during nights in which the 
bright-star RMS value was smaller than 0.04 mag, while the rest of the data 
is marked in red points. The histograms of airmass and color terms are shown 
on the top and right of the figure, respectively. The histograms are normalized 
by the sum of the number of data points. 

culated based on data from all nights in which the RMS of the 
best-fit residuals of bright stars is smaller than 0.04 mag and 
the limiting magnitud^is fainter than 20.5. We note that the 
/?-band color terms are roughly the same for all CCDs, while 
the g-band color terms have a wider range. 

8. COLOR TRANSFORMATIONS 

Approximate relations between the non-color-corrected 
PTF magnitude system and the SDSS magnitude system are 
given in Equations |4}{5] The exact relations between the two 
systems depend on the details of the object's spectral energy 
distribution and the atmospheric conditions at the time of the 
observations. 

The limiting magnitude is estimated based on the calibrated photometric 
parameters, the sky background and assuming a 3-(T detection limit. This 
property is stored in the image headers, as well as in the PTF database. 
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Fig. 10. — Seasonal variations in Oa/; and Oc,/; as a function of month. This 
plot is based on data taken between May 2010 and May 2011. 

Using Equation|4]and the color transformations between the 
Johnson-Cousins {UBVRJc) system and the SDSS systenQ 
we find an approximate transformation between PTF R-h&nd 
magnitude and the Johnson-Cousins system: 

R^F^R, + (/?c -/c) (0.3088 - 1.0517ac,«) 
-1-0.0704 -f0.2504ac« 

= Tij{{R,-h) + T2,R, (6) 

and similarly using Equation|5]we get for PTF ^-band: 

gPTF = y + (y - ^c) (0.9556 - 1 .652 1 a,,g) 
-0.0853-0. 1541 Ocg 

= ri,j(y-7?e) + r2,g. ' (?) 

Here T\ f and T2.f are the color-transformation coefficients for 
filter /. These parameters are listed in Table [3] 

Assuming] V = +0.03 and V - Rc = Rc - h =0,wt find 
for the AOV star Vega: 



pVega , 
"pTF 



:0.17-F0.23ac.ff, 



(8) 



http://www.sdss.org/dr6/algorithms/sdssUBVRITransform.html. 
http://www.sdss.org/dr6/algorithms/sdssUBVRITransform.html. 
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TABLE 3 

Median photometric-calibration parameters and associated parameter dispersion 
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0.008 


0.006 
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02 


0.202 


0.008 


-0.112 


0.025 


0.007 


0.006 


0.096 


0.121 


04 


0.219 


0.010 


-0.119 


0.027 


0.003 


0.006 
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05 


0.227 


0.009 


-0.117 


0.027 


0.002 


0.007 


0.070 


0.127 


06 


0.225 


0.008 


-0.106 


0.026 


0.005 
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0.121 
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11 


0.218 
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0.025 
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0.007 


0.080 


0.125 


All 


0.212 


0.014 


-0.112 


0.025 


0.005 


0.007 


0.086 


0.124 



Note. — The first column indicates the CCDID, where "AH" represents all the CCDs. Parameters Ti and T2 for band / are defined in Equations [6ll7l A represents the 68th percentile 
range divided by 2, a robust measure of scatter in a parameter. 



and 



^Pxf- -0.08+ 0.25 Ocg. 



(9) 



For the absolute magnitude of the Sun, assuming V = 
+4.82, V -Rc = +0.36, Rc-Ic = +0.32 and using SDSS 
color transformations, we find: 



?s 4.68-0. 11 Ocff, 



and 



g|^?.^5.12-0.44ac,, 



9. CONCLUSION 



(10) 
(11) 



A wide variety of scientific applications require photomet- 
ric calibration of astronomical images. In this paper, we pro- 
vide a description of the current PTF photometric calibra- 
tion process and products. We show that roughly half of the 
PTF data were taken under conditions that allow calibration 
of the images to an accuracy of 0.02-0.04 mag relative to 
SDSS. Most importantly, the rest of the data, taken under non- 
photometric conditions, can be calibrated using data taken in 
photometric conditions. At the bright end (« 15 mag), the 
PTF photometric calibration is stable to a level of ss 0.02 mag. 
Future improvements to the PTF photometric calibration pro- 
cess that are important for the PTF type-la and Hp supernovae 



cosmology projects will be described elsewhere. 

The PTF photometric calibration allows us to measure the 
calibrated magnitude of transients outside the SDSS footprint. 
We have used the PTF dataset to compile a photometric cat- 
alog. The first version of this catalog is described in Ofek et 
al. (in prep.). We note that our tests show that PTF data can 
deliver relative photometry, in single images, with precision 
as good as 3mmag (e.g., van Eyken et al. 2011; Agiieros et 
al. 201 1). This capability is being built into the PTF relative 
photometry pipeline and will be described in Levitan et al. (in 
prep.). 
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